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I  INTRODUCTION 


That  cosmic-ray  ionization  is  a  source  for  NO  in  the  stratosphere 

x 

,  * 

was  noted  by  Warneck.  He  and,  to  a  more  extensive  degree,  Brasseur 

and  Nicolet,2  examined  the  worldwide  contribution  of  this  NO  source 
’  x 

compared  with  other  sources,  and  it  appears  to  be  relatively  minor. 

Her.ce,  the  catalytic  destruction  o'  ozone  h  cosmic-ray-induced  NO^  pre¬ 
sumably  plays  a  small  role  in  dote  -mining  stratospheric  ozone  abundances. 
But  because  of  varying  solar  magnetic-field  effects  dui  ing  the  (~ll-yr) 
sunspot  cycles,  low-energy  cosmic  rays  incident  upon  the  higher  (^60  ) 
latitude  regions  of  the  earth's  atmosphere  are  modulated  with  the  sunspot 
(quasi)  period.  The  11-year  variation  in  stratospheric  ionization  has 
been  measured  for  many  years.  The  ionization  produces  N  and  N  ;  these 
in  turn  react  to  produce  predict. \ble  amounts  of  NO  .  The  abundances  and 
temperatures  in  the  stratosphere,  together  with  measured  reaction  rates, 
iiVly  that  essentially  all  of  these  free  N  atoms  and  ions  ultimately  re¬ 
sult  in  the  formation  of  NO.  Thus  there  is  an  oscillating  source  of  NO 
of  known  strength  and  spatial  ..istribution  in  the  stratosphere  at  the 
high  latitudes  of  both  hemispheres.  If  ..0  has  the  expected  qualitative 
effect  on  ozone  abundance,  some  cyclic  variation  in  stratospheric  ozone 
should  exist  with  essentially  the  periodicity  of  the  sunspot  cycle  (or 
more  precisely  that  of  the  modulation  of  low-energy  cosmic-ray  flux 
associated  with  it).  We  show  below  in  Section  IV  that,  at  those  two 
ozone  measuring  stations  (Tromso  and  Arosa)  where  published  records  ex¬ 
tend  back  40  years,  the  measured,  integrated  ozone-column  densities  not 


★ 
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only  appear  to  have  the  expected  periodicity  but  also  show  the  calcu¬ 
lated  phase  lag  relative  to  the  sunspot  cycle.  Worldwide  data,  available 
after  tne  mid-1950's,  show  a  dependence  of  this  correlation  on  latitude, 
and  the  time  lag  and  relative  amplitude  nt  different  stations  are  in 
reasonable  agreement  with  the  calculated  ones.  The  absolute  magnitude 
of  the  periodic  ozone  modulation  is  also  shown  to  be  consistent  with  that 
expected  from  typical  models  for  the  stratosphere,  wherein  NO  plays  a 
major  role  in  controlling  the  ozone  abundance.  These  various  coincidences 
between  the  observ  i  and  calculated  quasi-periodic  modulations  of  ozone 
abundance  suggest  that  the  known  variation  of  stratospheric  NO  injection 
produces  the  measured  variations  of  ozone.  Having  established  that  the 
ozone  variability  is  roughly  proportional  (with  some  time  delay)  to  the 
change  in  NO  density,  then  the  effect  of  known  NO  injection  into  the  same 
regions  of  the  stratosphere  by  SST  aircrait  (or  nuclear  explosions)  is 
calculable  (Section  7). 


II  SUNSPOT-CYCLE  MODULATION  OF  STRATOSPHERIC 
COSMIC-RAY  IONIZATION 

The  formation  of  ion  pairs  by  cosmic  rays  has  been  studied  at  high 
latitudes  as  a  function  of  altitude  over  three  sunspot  cycles. 4  Re¬ 
sults  are  summarized  in  Figures  1  through  4.  It  is  clear  that  ionization 
deposited  at  stratospheric  heights  is  strongest  over  the  polar  caps, 
where  the  magnetic  field  is  weakest,  and  that  a  sharp  drop  occurs  below 
60°  geomagnetic  latitude.  The  pe^k  ion  production  is  between  11  and  15 
km  above  the  Arctic  tropopause  at  8  to  9  km.  The  mean  peak  is  about  35 
ion-pair  cm  s,  wit'  an  .mplitude  through  the  sunspot  cycle  of  about 

7 

±15  percent.  The  total  ionization  in  a  column  averages  about  4.5  X  10 
-3  -1 

ion-pair  cm  s  .  with  an  amplitude  through  the  sunspot  cycle  of  about 

7 

±li>  percent.  The  to+al  ionization  in  a  column  averages  about  4.5  X  10 
-2  -1 

ion -pair  cm  s 
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IONS  cm*3  T1  ATM*’  OF  AIR  ZU"'CH  SUNSPOT  NUMBER 


FIGURE 


1  CORRELATION  OF  SUNSPOT  NUMBER  WITH  IONIZATION  AT 
SELECTED  PRESSURES  OVER  THULE,  GREENLAND.3  The 
pressures  are  specified  by  equivalent  altitude  in  a  Standard  Mode. 
Atmosphere.  (The  ionization  must  be  multiplied  by  the  pressure 
in  atmospheres  to  find  ion  pairs  crri"3.) 
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SOLAR  MINIMUM 
_ (1954  and  19651 


SOLAR  MAXIMUM 
(1958  and  1969) 


ION  PAIR  PRODUCTION  — cm"3  s"1 

FIGURE  2  VARIATION  OF  COSMiC-RAY-INDUCED  ATMOSPHERIC  IONIZATION  IN  THE 
ATMOSPHERE  OVER  THULE,  GREENLAND,  BETWEEN  SUNSPOT  MAXIMUM 
AND  MINIMUM.  Source:  Ref.  2,  from  data  of  Ref.  4. 


MEAN 

TROPOPAUSE 


II 

(ion  pairt  \ 
cm3.  / 


FIGURE  3  AMPLITUDE  OF  OSCILLATING  COMPONENT  OF  '.TRATOSPHERIC 
IONIZATION  FOR  THE  DATA  OF  FIGURE  1  AND  THE  STANDARD 
MODEL  ATMOSPHERE  AS  A  FUNCTION  OF  ALTITUDE  (km).  The 
abscissa  i*  ’he  maximum  ionization  variation  21 1;  the  ionization  I  = 

I0  +  Ii  cos  wt  in  ion  pairs  cm~3  y-1. 
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GEOMAGNETIC  LATiTUCT 


FIGURE  4 


IONIZATION  vs.  LATITUDE  FOR  CONSTANT 
OVERHEAD  AIR  MASS  OF  15  g  cm'2  FOR 
THREE  DIFFERENT  YEARS.  In  1937  the 
sun  was  at  a  maximum  of  activity  while  in 
1954  it  was  at  a  minimum.  The  pressure  of 
15  g  cm"2  corresponds  to  about  30  km 
altitude.  The  vertical  arrows  show  the 
latitudes,  at  two  zenith  angles,  for  cutoff 
by  the  geomagnetic  field  of  primary  cosmic 
rays  that  can  just  penetrate  an  air  mass  of 
15  g  cm'2.  Source:  Ref.  5. 


Ill  MODULATED  NO  INJECTION  BY  COSMIC  RAYS 


Calculations  by  Dalgarno6  on  the  distribution  of  collision  products 
iron  ionizing  particles  in  air  were  used  by  Warneck1  to  estimate  that 
0.33  NO  molecule  is  produced  for  each  ion  pair.  However,  Brasseur  and 
Nicolet2  argue  that  the  dissociation  of  1'  by  secondary-electron  impact 
should  be  increased  by  a  factor  of  three,  and  the  NO  productir  1  should 
be  raised  by  about  the  same  factor.  More  detailed  calculations  by 
Gilmore^  give  a  yield  of  1.3  N  atoms  per  ion  pair. 

This  is  also  effectively  the  rate  of  production  of  NO  per  ion  pair 
produced  by  cosmic  rays  via  reaction  (2)  below.  In  a  note  added  in  proof 
Warneck  mentions  that  N  can  also  destroy  NO.  Indeed,  the  reaction 

N  +  NO  -  N2  +  0  (1) 

has  one  of  the  fastest  two-body  rate  coefficients  known, 8,9  The  rate  has 

-11  3  -1 

been  measured  as  about  2  X  10  cm  s  at  300  K  and  it  varies  as 
5  X  10  11  exp  (-100/T)  at  somewhat  higher  temperatures  (475°  to  750  K) . 
Hence  we  adopt 


-11  3  -1 

k^  =  2.8  X  10  exp  (-100/T)  cm  s 

for  stratospheric  temperatures. 

The  NO  creation  mechanism  is 

N  +  0  -  NO  +  0  (2) 

2 

7 


4 


Vlastraras  and  Winkler,10  and  Wilson11  give  formulae  yieldirg  results 
that  differ  by  a  factor  of  two  when  extrapolated  to  the  200  K  regime. 

The  Vlastraras-Winkler  coefficient,  which  gives  the  larger  values,  is 

•12  3  -1 

k  =  6.5  X  10  exp  (-3500/T)  cm  s 
2 

Figure  5  gives  [NO]  fractional  concentrations  for  three  different  formulae 

for  the  rate  coefficient,  k^,  when  N  and  NO  are  in  chemical  equilibrium 

through  reactions  (1)  and  <2).  In  the  lower  stratosphere,  where  the 

—  18  —  3 

pressure  is  0.1  atm,  the  density  is  [M]  =  ?  X  10  cm  and  T  =  220  K, 

-9 

the  equilibrium  concentration  [N0]/[M],  exceeds  4  x  10  .  The  equilibrium 

mole  fractions  of  NO  in  the  presence  of  a  continued  injection  of  N  atoms 

are  plotted  versus  altitude  in  the  Standard  Model  Atmosphere*2  in  Figure 

6.  The  actual  relative  concentration  of  stratospheric  NO  in  a  steady 

state  is  not  well  known.  While  theoretical  models  suggest  chemical- 

-9 

equilibrium  concentrations  exceeding  10  ,  measurements  around  21  km  by 

Ridley  et  al.la  give  only  about  10  10.  However,  spectroscopic  measure- 

-9 

ments  by  Totl  et  al.14  from  aircraft  indicate  concentrations  of  10  in 

the  altitude  range  11  to  26  km.  Even  if  the  equilibrium  NO  abundance 

for  thermalized  chemistry  were  achieved,  N  atoms  produced  from  cosmic-ray 

ionization  would  still  make  NO  more  rapidly  than  it  would  destroy  it. 

Gilmore7  finds  that  more  than  half  of  the  N  atoms  are  created  in  the 
2  2 

metastable  states  D(2.37  eV)  or  P(3.56  eV),  which  are  not  readily 

quenched  by  thermal  collisions  with  N  molecules.  With  this  energy 

z 

available  to  counteract  the  large  activation  energy  of  reaction  (2), 

5 

the  rate  coefficient  k  could  be  increased16’ 16  by  a  factor  of  10  .  (In 

addition,  the  initial  kinetic  energy  of  the  N  atoms  will  likely  exceed 

the  local  thermal  energies  and  remain  so  for  the  first  few  collisions 

with  N  molecules.)  Thus,  more  than  half  the  N  atoms  are  expected  to 
2 
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200  220  240  260  280 

TEMPERATURE  — °K 

FIGURE  5  PARTIAL  MOLE  FRACTION  OF  NO  IN  PARTS  PER  BILLION  (10-9)  VERSUS 
TEMPERATURE  AT  WHICH  THE  NET  RATE  FOR  PRODUCTION  OF  NO  BY  N 
VANISHES  —  i.e.,  N  +  NO  -*  N2  +  0  =  N  +  02  -*•  NO  +  0,  WHERE  N  IS 
ASSUMED  TO  BE  THERMALIZED  IN  THE  GROUND  STATE.  The  three  curves 
represent  three  extrapolations  from  measurements  at  higher  temperatures  (cf.  Section  III). 
k2a  *  6.5  x  10-12  exp(-3500/T),  source:  Ref.  10;  k2b  *  2.33  x  10-1 1  exp(-39i»0/T), 
source:  Ref.  11;  k2c  *  4.98  x  10-13  j'A  exp(-3565/T),  source:  Ref.  11. 


HEIGHT  —  km 


FIGURE  6 


PARTIAL  MOLE  FRACTION  OF  NO  IN  FARTS  PER  BILLION  FOR  WHICH 
THERMALIZED  GROUND  STATE  N  ATOMS  WOULD  PRODUCE  NO  NET 
NO  IN  A  STANDARD  MODEL  ATMOSPHERE.  The  maximum  and  minimum 
rate  coefficients  of  Figure  5  are  used. 


make  NO  through  reaction  (2)  no  matter  how  high  the  NO  fractional  abun¬ 
dance.  In  our  numerical  estimates  we  shall  assume  all  do. 

Throughout  the  above  discussion  we  have  supposed  for  definiteness 

that  NO  is  entirely  in  the  form  NO.  However,  if  it  were  largely  NO  , 
x  2 

as  it  may  be  in  the  lower  stratosphere,  our  conclusion  that  N  mainly 
anufactures  NO^,  instead  of  destroying  it,  is  unchanged.  The  rate  co¬ 
efficient  for  N  +  NO  reactions  at  300°  K  is  also  large,17  about 
o  -11  3  -1  2 

1.8  x  10  cm  s  However,  there  are  four  possible  end  products  and 

the  creation  of  2  NO's  has  a  1/3  branching  ratio.  Hence  the  destruction 

of  NO  ic  down  by  a  factor  of  about  3,  when  the  reactant  is  NO  ,  compared 
x  2 

with  the  loss  rate  of  NO  by  reaction  (1).  Consequently,  any  NO^  present 

will  have  the  net  effect  of  increasing  the  NO  caused  by  cosmic-ray- 

x 

produced  N. 
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IV  EVIDENCE  FOR  THE  CORRELATION  OF  OZONE 
WITH  THE  SOLAR  CYCLE 

Claims  for  the  existence  of  a  correlation  of  ozone  abundance  with 
sunspot  activity  have  an  extensive  history.  Recent  analyses  of  30-month 
running  averages  of  ozone  column  densities  by  Angell  and  Korshover18 
seem  to  offer  significant  support  for  such  a  correlation.  Their  processed 
data  are  given  in  Figures  7  and  8.  Only  two  stations,  Tromso  and  Arosa, 
are  given  with  data  extending  back  far  enough  to  cover  several  sunspot 
cycles.  Both  stations  appear  to  have  very  significant  correlation  with 
solar  magnetic  activity.  For  the  sunspot  cycle  around  1960,  after  the 
extensive  worldwide  ozone  data  collection  of  the  mid- 1950's  began,  hemi¬ 
spheric  and  equatorial  ozone  also  seem  to  correlate  with  sunspots  but 
with  smaller  amplitude.  (The  30-month  running  average  Northern  Hemisphere 
data  seem  to  be  irregular  in  the  interval  around  1960;  this  appears  not 
only  in  the  extensive  data  of  Figure  8  but  apparently  also  in  the  more 
erratic  ozone  data  from  Arosa  and  Tromso.  An  effect  on  the  30-month 
averaging  of  the  extensive  Arctic  nuclear-explosion  program  beginning 
mid-1961  might  be  expected  from  some  ozone  suppression  by  bomb-injected 
NO.  If  so,  the  data  suggest  the  averaged  effect  may  be  of  the  order  1/2 
percent  for  the  Northern  Hemisphere.  During  this  particular  sunspot 
cycle,  Northern  Hemisphere  correlation  data  may  be  suspect.)  If  the 
correlation  exists,  the  effect  is  strongest  near  the  poles  (Tromso,  70°) 
where  the  ozone  maximum  typically  lags  the  sunspot  maximum  by  two  to 
three  years;  it  is  progressively  weaker  and  has  a  greater  time  lag  at 
lower  latitudes.  Just  this  qualitative  behavior  would  be  expected  from 
cyclical  injection  at  high  latitudes  of  some  catalyst  which  destroys 
ozone  (or  enhances  it — depending  on  the  phase  relation  of  catalyst 


12 


TOTAL  OZONE  VARIATION  —  pwcant 


TROMSO 


AROSA 


38  MONTHS 


TROPICS 


1960  1966 


FIGURE  7  COMPARISON  OF  TEMPORAL  VARIATION  IN  SUNSPOT  NUMBER  (top)  WITH 
THE  PERCENTAGE  DEVIATION  FROM  THE  MEAN  OF  THE  30-MONTH 
RUNNING  AVERAGES  OF  TOTAL  OZONE.  The  data  for  the  tropics  are 
"smoothed"  and  have  a  scale  one-tenth  that  of  Arosa  (47°N)  and  Tromso  (70°N). 
For  Tromso  and  Arosa,  for  which  over  30  years  data  are  available,  the  average 
phase  lag  of  ozone  variation  relative  to  *hat  of  sunspot  count  is  given  in  months. 
Source:  Ref.  18. 
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FIGURE  8  COMPARISON  OF  THE  TEMPORAL  VARIATION  IN  SUNSPOT  NUMBER  WITH 
THE  SMOOTHED  PERCENTAGE  DEVIATION  FROM  THE  LINEAR  TREND  OF 
30-MONTH  RUNNING  AVERAGES  OF  TOTAL  OZONE.  Fits  for  the  lag  in 
response  of  the  world  zone  and  the  Northern  Hemisphere  zone  relative  to  the 
sunspot  variation  are  given  at  right.  Source:  Hef.  18. 
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injection  with  sunspot  cycle):  the  relative  time  lag  and  the  amplitude 
decrease  with  increasing  distance  from  the  poles  are  characteristic  of 
eddy  diffusion  away  from  an  oscillating  polar  source.  In  the  next  sec¬ 
tions  we  show  how  the  observations  are  more  quantitatively  consistent 
with  the  model  calculations.  If  the  solar  cycle  were  to  modulate  strato¬ 
spheric  ozone  by  altering  solar  radiation  in  an  appropriate  part  of  the 
ultraviolet  spectrum,  its  initial  effect  on  the  earth  would  probably  be 
worldwide  and  simultaneous,  thereby  giving  a  temporal  and  spatial  pattern 
to  resulting  ozone  variations  very  different  from  that  found  by  Angell 
and  Korshover  18 
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V  CALCULATED  TIME  DELAYS  AND  RELATIVE  AMPLITUDES 
MODUI  VITD  OZONE  AT  VARIOUS  LATITUDES 


We  consider  an  idealisation  of  the  stratosphere  in  which  the  resi¬ 


dence  time  t  for  NO  is  (cf.  Figure  9) 


t  ~  2  to  6  years 


and  in  which  horizontal  diffusion  is  described  by  r  'onventional  constant 


horizontal  eddy-diffusion  coefficient, 
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FIGURE  9  "MEASURED"  STRATOSPHERIC  BURDEN  OF  BOMB  MANUFACTURED  C14. 

Source:  Refs.  19  and  20.  The  left  scale  is  that  for  C^4.  The  right  is  for 
NO  on  the  assumption  that  a  one  megaton  air  explosion  which  locally  forms 
2  x  1026  c’4  atoms  (Ref.  21)  also  forms  <  x  1031  NO  atoms.  The  latter 
is  the  lower  bound  of  Gilmore  (Ref.  7)  and  almost  that  of  Ref  21.  The 
solid  line  is  the  total  stratospheric  burden;  the  dashed  line  is  the  30-month 
running  average  for  the  Northern  Hemisphere  alone.  Bomb  yield,  are  given 
in  megatons  for  testing  intervals. 


Then  in  the  presence  of  an  NO  source  J,  we  have 


J 


[NO  j 

T 


2 

KV  [NO] 

Cd 


(5) 


where  V  is  the  two-dimensional  Laplacian  over  the  surface  of  the  ter- 
2 

restrial  sphere.  We  assume  the  source  has  an  oscillating  component 


■J  =  J 


o 


iiut, 

4-  J  e 
1 


(6) 


with 


in 


2tt 

11  years 


(7) 


T’tv?  amplitude  is  constant  in  the  high  latitude  regions,  Q  &  ~  60  , 

and  is  zero  away  from  the  polar  regions  where  9  <  <5^.  Equations  (5)  and 

(6)  are  easily  solved  in  terms  of  the  Legendre  functions  (cos  Q)  and 

Q  (cos  Q)  but  with  complex  v.  The  analogous  solution  on  an  infinite 
v 

plane  with  uniform  oscillating  injection  within  a  circle  of  radius  r^  is 
developed  in  Appendix  A.  For  the  special  case 

.1/2  r2 

(u>  +  T  )  ^  «  1  (8) 


this  solution  becomes  [Eq.  (A-20)j 

2 

0,4  J^r^ 

[NO] - - - -  exp  [ itu(  t  -  A>]  (9) 

K 

for  r<  r  with  the  time  lag  A.  For  the  parameters  of  Eq~  (3),  (4), 
o 

and  (7),  we  obtain 

A  ~  1.3  years  .  (10) 
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t  (Tiomso)  ~  2.3  years 
L 


(15) 


We  assume  further  that  total  or, one  abundance  responds  to  altered  NO  abun¬ 
dance  in  a  time  much  less  than  the  various  time  lags  between  sunspot 
maxima  and  NO  minima.  From  Figure  7,  a  best  fit  of  Angell  and  Korshover's] 
compilation  of  ozone  data  gives  1  (Tromso)  ~  2.7  years,  in  satisfactory 

i-l 

agreement  with  the  model  calculation.  For  Arosa  (0  *  47°)  theory  gives 

t  (Arosa)  -  tT  (Tromso)  ~  0.3  years  ,  (16) 

compared  with  the  published  observational  fit  of  38  -  32  months  =  0.5 
years.  The  calculated  amplitude  at  Arosa  relative  to  Tromso  of  0.6  is 
also  close  to  the  "observed"  one.  For  equatorial  regions  the  model 
suggests 


t  (equatorial)  -  t,.  (Tromso)  ~  1.1  years  ,  (17) 

compared  with  an  observed  lag  of  rbout  2.5  years.  The  computed  amplitude 

in  the  tropics  relative  to  that  at  high  latitudes  is  about  0.2,  perhaps 

about  twice  that  observed.  (Better  fits  to  the  time  lags  are  obtained 

with  a  somewhat  smaller  horizontal  diffusion  coefficient  K;  the  relative 

-1/2 

time  lags  for  0<  0  are  proportional  to  K  .)  The  rough  agreement 

o 

with  the  Angell-Korshover  analyses  gives  some  additional  support  to  the 
reality  of  the  observed  sunspot-ozone  correlation  v.id  to  cosmic-ray  modu¬ 
lated  NO  in  the  high  latitude  regions  as  its  origin.  The  quantitative 
agreement  with  relative  amplitudes  is  somewhat  fortuitous  because  of  the 
varying  NO  altitude  distributions  at  different  latitudes  (cf.  Section  VII) 


VI  ESTIMATES  FOR  ABSOLUTE  MAGNITUDE  OF  COSMIC- RAY- INDUCED 

MODULATIONS  OF  OZONE 


Over  the  polar  regions  (0  >  60°)  the  calculated  [NO]  is  given  by 

Eq.  (9)  with  A  now  replaced  by  t  of  Eq.  (13).  Then  for  the  change  in 

L 

the  column  density,  J  [NO]  dz,  we  have 

A  J  [no]  dz  ~  1  X  10  cos  [uJ(t  -  2.5  years)]  , 

for  0  >  60°  .  (18) 


From  the  magnitude  and  altitude  dependence  of  cosmic-ray-modulated  ioniza 
tion  of  Figure  3,  we  obtain  a  column  production  of 


J  ~  8 


X  10  molecules  cm  s 


(19) 


and 


A  J  [NO]  dz  ~  8  X  10  cos  uu(t  -  2.5  years)  molecules  cm  , 

(0  >  C0°)  .  (20) 

The  corresponding  ozone  variation  measured  at  Tromso  is  about  ±5  percent 
(cf..  Figure  7). 

A  similar  estimate  follows  for  worldwide  averages.  The  amplitude 
of  the  worldwide  average  of  modulated  NO,  A  j*  [NO]  dz,  is  independent  of 
the  horizontal  diffusion;  eliminating  the  K  term  from  Eq.  (5),  we  find 


) 


(21) 


L  J  [NO]  dz 

with  a  mean  phase  lag  of 


J,  cos  ( cut  -  i)  ur 
l _  o 


+  T 


-2 


2ttR 


e 


4  =  tan  UOT 


Then 


[NO]  dz  ~  4  X  1013 


cos  (UDT  -4  ) 


(22) 


(23) 


for  t  ~  3  years. 

The  corresponding  change  in  amplitude  for  the  cyclic  part  of  the 
world  or  Southern  Hemisphere  ozone-column  density  from  Figure  8  is  about 
0.5  percent,  consistent  with  a  bigger  effect  from  about  a  fivefold  larger 
NO  perturbation  at  Tromso  (with  a  different  altitude  dir tribution) . 
Appendix  B  summarizes  arguments  for  a  very  rough  dependence  of 


[°3] 


3  AJNOJ 
8  [NO] 


(24) 


in  stratospheric  models  (e.g. ,  Johnston,23  and  Chang  and  Duewer24),  where 
ambient  [NO]  reduces  ozone  to  about  half  what  it  would  be  from  the  Cnapman 
cycle  alone.  Then,  from  Eqs.  (20)  and  (23)  and  the  corresponding  observed 
ozone  oscillations,  the  ambient  NO  column  density  would  be 


J  [NO]  dz  ~  |  y 


8  x  10 
5  X  10' 


13 


=  6  X  10 


14  molecule 


(25) 


cm 


9  _3 

corresponding  to  about  10  NO^_  molecules  cm  in  the  lower  stratosphere 
or  a  relative  concentration  of  about  1  ppb.  This  is  lower  but  roughly 
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of  the  same  order  of  magnitude  as  assumed  in  models  in  which  NO  contributes 
significantly  to  the  control  of  ozone  (cf.  Figure  10,  Chang  and  Duewer,24 
and  Johnston  et  al.20).  We  emphasize  that  the  estimate  of  Eq.  (25)  is 
extremely  rough  and  can  be  taken  only  as  support  for  the  proposition 
that  the  computed  cosmic-ray-modulated  change  in  ozone  has  a  reasonable 
magnitude. 

Thus  the  magnitude  as  well  as  the  time-lag  and  latitude-dependence 
of  ozone  variations  seem  consistent  with  the  model  for  cosmic-ray- 
modulated  NO  as  the  origin  of  the  11-year  cyclic  ozone  variation. 


FIGURE  10  AN  OZONE  PROFILE  OBSERVED  AT  PANAMA 
(9°N)  COMPARED  TO  CALCULATED  ONES  FOR 
THE  EQUATOR  WITH  NOx  MOLE  FRACTIONS 
OF  3.  10,  AND  100  ppb  TO  ACCOUNT  FOR  THE 
OZONE  DEFICIT.  Source:  Ref.  23. 
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VII  OZONE  DEPRESSION  FROM  STRATOSPHERIC  NO  INJECTION 
BY  ARTIFICIAL  SOURCES 


If  the  relation  between  sunspot  activity  and  ozone  suppression  is 
real  and  caused  by  the  calculated  varying  stratospheric  NO  injection, 

+hen  existing  ozone  data  can  be  used  to  predict  suppression  from  other 
sources  of  NO. 

A.  Nuclear  Explosions 

Conventional  analyses  of  expected  NO  stratospheric  injection  from 

34 

nuclear  explosions  give  1  to  3  x  10  NO  molecules  injected  into  the 
stratosphere  during  the  14-month  period  beginning  August  1961  (Foley  and 
Ruderman,22  and  Gilmore7).  About  90  percent  of  the  explosion  yields 
were  deposited  at  Novaya  Zemla  near  the  North  Pole  (0  =  72°).  The  alti¬ 
tude  distribution  of  this  NO  is  presumably  very  similar  to  that  of  bomb- 
14 

created  C  .  Some  distributions  inferred  from  sampling  data  are  given 

in  Figures  11  through  16  (Telegdas19 ) .  It  appears  that  C14  in  the  polar- 

cap  region  (0  >  60  )  has  an  initial  distribution  which  peaks  at  around 

18  km  aid  a  full  width  of  about  10  km.  This  is  qualitatively  not  very 

different  from  the  computed  altitude  distribution  of  the  modulated  cosmic- 

ray-injected  NO  of  Figure  3;  therefore,  the  bomb-injected  NO  (assuming 

14 

it  follows  bomb-injected  C  )  may  be  expected  to  spread  with  an  altitude 
distribution  at  each  latitude  similar  to  the  distribution  of  NO  from  the 
11-year  oscillating  polar  source.  That  is,  a  direct  comparison  between 
the  magnitudes  of  the  ozone  effects  from  the  two  different  kinds  of  polar 
sources  should  be  feasible. 
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FIGURE  14  RESULTS  OF  ATMOSPHERIC  C14  SAMPLING  AND  EXTRAPOLATION  — 
MARCH-MAY  1963  (Telegdas1 9 ) .  (The  measured  numbers  and  the  isodensity 
line  labels  are  in  units  of  105  C14  atoms  per  gram  of  air.) 
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FIGURE  '16  RESULTS  OF  ATMOSPHERIC  SAMPLING  AND  EXTRAPOLATION  — 
MARCH-MAY  1966  (Telegdas*  *).  (The  measured  numbers  and  the  isodensity 
line  labels  are  in  units  of  10^  C^  atoms  per  gram  of  air.) 


Johnston,  Whitten,  and  Birks20  have  normalized  bomb- injected  NO  to 

90 

observed  bomb-manuf ac tured  Sr  .  A  similar  analysis  based  upon  bomb- 
14. 

manufactured  C  is  given  in  Figure  9.  A  one  megaton  explosion  produces 
16  14  4 

2  x  10  C  nuclei  (within  about  3  x  10  cm  of  a  low-altitude  explosion). 

According  to  Gilmore  the  minimum  NO  production  from  the  same  explosion  is 
31 

about  4  X  10  molecules.  The  radii  of  the  respective  volumes  within 

which  each  species  is  made  is  comparable.  Thus  a  reasonably  conservative 

5  14 

estimate  suggests  2  x  10  NO  molecules  per  C  atom.  Figure  9  gives  the 

14  X 

total  measured  excess  C  burden  of  the  stratosphere  from  bomb  explosions. 

28  14  33 

This  peaks  at  3.6  X  10  C  ,  implying  7  X  10  NO  in  the  first  quarter 

of  1963.  The  30-month  running  average  for  the  Northern  Hemisphere  alone 

28  14 

stays  at  near  2  x  10  C  nuclei  from  mid-1962  to  the  end  of  1963, 

33 

corresponding  to  perhaps  4  x  10  NO  atoms  in  the  Northern  Hemisphere 

27 

stratosphere  during  this  interval.  Since  there  were  about  5  X  10 
14 

excess  C  atoms  in  the  Northern  Hemisphere  stratosphere  from  previous 

tests  before  the  intense  testing  of  1961-62,  the  extra  NO  inserted  might 

33 

conservatively  be  estimated  as  averaging  at  least  3  x  10  NO  during  1963, 

15  -2 

This  gives  an  average  column  density  of  1  X  10  cm  with  about  the  same 

altitude  distribution  as  modulated  NO  from  the  cosmic  rays.  (This  is 

15 

less  than  the  total  minimum  NO  of  4  x  10  obtained  by  averaging  total 

produced  NO  because  of  losses  from  the  stratosphere,  30-month  averaging 

and  perhaps  imperfect  injection.)  The  estimates  of  Section  VI  imply  that 

14  -2 

a  change  in  average  NO  column  density  of  10  cm  produces  a  1  to  2  per¬ 
cent  change  in  ozone.  The  lower  estimate  together  with  the  "minimum"  NO 
bomb  yield  then  suggest  that  30-month  running  averages  of  total  ozone 
should  show  a  drop  of  3  percent  by  mid-1961,  continuing  to  a  10  percent 
depression  by  mid-1962,  an  1  remaining  at  that  level  until  the  end  of 

1963.  From  analyses  of  worldwide  ozone  data,  Johnston,  Whitten,  and 

# 

Birks30  claim  evidence  for  a  world  average  ozone  decrease  of  3.3  percent 
during  the  1960-62  interval  of  extensive  testing.  Such  a  decrease  does 


30 


not  appear  in  the  analysis  of  Angell  and  Korshover18  who  find  "little 
evidence  of  a  reduction  in  total-ozone  due  to  nuclear  testing"  down  to 
less  than  a  percent. 

Therefore,  if  yields  of  the  bomb  NO  into  the  stratosphere  are  as 
large  as  those  estimated  by  Gilmore7  and  by  Foley  and  Ruderman,22  there 
seems  to  be  a  conflict  with  our  proposed  explanation  of  an  11-year  ozone 
cycle  and  possibly  also  with  ozone  observations  during  and  after  the 
1961-62  period  of  heavy  testing.  A  much  smaller  amount  of  bomb-injected 
stratospheric  NO^  in  a  form  in  which  it  catalyzes  ozone  destruction  would 
be  comforting  to  theories  for  NO  suppression  of  ozone. 


B.  Superset io  Transports 

If  fleets  of  supersonic  transports  were  to  fly  only  at  latitudes 
greater  than  60°  with  the  altitude  distribution  of  Figure  3,  then  a  com¬ 
parison  of  the  effects  from  the  sources  of  Eqs.  (20)  and  (23)  with  those 
of  aircraft  would  involve  no  other  additional  assumptions  except  that 
the  ozone  response  is  linear  with  NO.  The  exhaust  NO  injected  into  the 
stratosphere  would  be  expected  to  diffuse  over  the  globe  with  the  pat¬ 
terns  and  timing  of  Figures  11  through  16.  Since  t  ~  2  to  6  years  is 

somewhat  greater  than  the  diffusion  time  from  the  pole  to  the  equator 

2  2  -1 
(tt  R^/ieK  ~  1  year)  a  steady  source  I  (NO  molecules  s  ),  half  in  each 

polar  cap,  would  raise  stratospheric  NO  roughly  uniformly  by 


where  I  „ 
33 

per  year, 
predicted 
altitudes 


It  .  ..  _„14  /  t  \ 

~  4  X  10  I  '  I  - 

.  2  \2  years/  33 

4ttR 


NO  molecules 


(26) 


cm 


e 

33 

is  the  NO  input  from  the  SST  fleet  in  units  of  10  molecules 
A  comparison  with  the  effects  of  Eqs.  (20)  and  (23)  gives  a 
average  worldwide  suppression  from  an  SST  fleet  at  appropriate 
and  latitudes,  after  a  steady  state  is  achieved,  of 
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(4  to  8)  (  — -  )  I  percent,  if  the  effect  is  small  enough  so  that 

y  2  years  J  33 

a  linear  extrapolation  remains  valid.  Johnston25  has  estimated  that  500 
SST’s  flying  seven  hours  per  day  vould  give  I  ~  30,  so  that  the 
corresponding  effect  on  ozone  would  be  very  large.  The  SST  flying 
altitude,  near  20  km,  is  above  most  of  the  cosmic-ray-modulated  NO 
insertion  in  the  polar  regions.  However,  the  naturally  varying  source 
is  well  above  the  low  polar  tropopause,  so  postulating  similar  results  at 
other  latitudes,  where  much  higher  tropopauses  ^ ill  have  NO  sources  just 
above  them,  may  be  reasonable  in  order  to  predict  the  effects  of  NO  from 
SST  .’lights. 
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Appendix  A 


NO  HEMISPHERICAL  VARIATION  FROM  OSCILLATING  POLAR  SOURCES 


itut 

In  a  plane  with  a  spatially  constant  source  J^e  within  a  circle 

of  r  <  r  and  zero  source  for  r  >  r  ,  the  solution  of  Eq.  (5)  is 
o  o 


[NO]  =  e 


iiut 


iiu  +  - 
T 


T  +  BI  (z) 
1  o 


( A-l) 


for  r  <  r  ,  and 
o 


[NO]  =  e 


iiut 


K<z)] 


(A- 2) 


for  r  >  r  .  In  these  equations 
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■  tt  (“" + ;) 


(A- 3) 


A  = 


K'(z  \l  (z 

_  _  °\  °/  °\  o 

(lu)  +  ”)  K0  (z0)  L  K°(Zo)Io(Zo 


-.-1 


(A- 4) 


and 


B  =  A 


(A-5) 


Here  I  and  K  are  the  usual  modified  Bessel  functions  and  z  is  the 
o  o  o 

value  of  z  defined  in  Eq.  (A-3)  at  the  source  radius  r  =  r  .  We  asso- 

o 

elate  this  radius  with  that  of  the  polar  cap  for  NO  Injection  measured 
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■ 


along  a  meridian  from  the  pole  to  latitude  60  .  Then 


8 

r  ~  3,3  x  10  cm 
o 


Then  with  t  ~  3  years,  we  have 


(  A—  ( 


(A-: 


Within  the  cap  and  at  z  =  z  we  therefore  keep  only  the  first  term  in 

o 

the  expansion 


I 

o 


(  A—  1 


Then  the  solutions  (A-i)  through  (A-5)  reduce  to 

io)t 
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o 


|iu)  +  (1  -  a)  Ko(Zo) 
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z  <  z 


z  >  z 


Further  approximations,  valid  in  the  limit  z  ->  0,  are 


(A-! 


(A-li 


(  A-l 
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( A-14) 


and 


Ko(zo)  ~  n  -  y  +  In  2 
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where  y  is  the  Euler  constar  0.57721  .... 


From  Eqs.  (A  11) 
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When  Eq.  (A-3)  for  z  is  substituted  into  Eq.  (A-17) 
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[NO]  ~ 


icut  2 
J  e  r 
1 _ o 

2K 


0.11  -  0.5  In 


f  Hi) 


z  <  z  «  1 
o 


With  valuer  for  m,  t  (=  3  years),  K,  and  r  obtained  from  Eqs 

o 

(7),  and  (11)  of  the  main  text,  we  find 
Jiro 

[NO]  ~  0.4  ——  cos  u>  ( t  —  1.3  years)  =  [NO]  ,  z  <  z 

K  o  o 


For  z  >  we  use  only  the  first  term  in  the  expansion 


(A-15) 

through 

(A-16) 

(A-17) 

( A-18) 

(A-19) 

.  (4),  (3), 

(A- 20) 
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1 


9 


•  • 


( A-21 ) 


K  (z) 
o 


•(=) 


.1/2 


+ 


82  2<8z)2 


Then  Eq.  (A-18)  becomes 
.  1/2 


[NO]  =  [NO] 


=  [NO] 


M  /K) 


ft 
.1/2 

(*o\  /  \  -1/2  /  2  -2X1  4 

exp[-(?  -  v  r  +T  ) 


16/2 


r  >  r 


and 


6  =  tan  urr  ~  1.0  rad 


for  t  ~  3  years  and  cp  the  radian  angular  distance  froji  the  pole, 
cal  substitution  gives 


[NO] 


cpQ ) ( 1 . 1  +  C. 62  i) 


( A-22) 


( A-23) 

Numeri- 

( A-24) 
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Appendix  B 


OZONE  VARIATION  FROM  CHANGE  OF  NO 


In  Section  VI  we  relate  the  variation  of  [NO]  to  that  in  [O  ] .  We 

3 

can  obtain  a  rough  approximation  for  that  relationship  without  a  detailed 
knowledge  of  reaction  rates,  except  to  invoke  the  assumption  that  NO  is 
necessary  to  account  for  lowering  [0  ]  by  a  factor  a  below  the  ambient 
abundance  indicated  by  pure-oxygen  photochemistry. 

Odd  oxygen  is  governed  by 

”  ([O]  +  [03]j  =  J  -  K^NO]  [0]  -  K2[0]  [03]  ,  (B-l) 


where  J  is  the  production  rate  of  odd  oxygen  by  photodissociation  and 


(B-2) 


These  relations  yield 


3  /  Vl/2 

(l  +  K5[NO]/£03]) 

The  factor  a  defined  above  means  that,  for  the  steady-state  abundances 
of  [NO]  and  [0  ]  , 

A  Q  A7 


[NO] 

Ks  Ttf  -  *  -  1  5  P  *  <B'4) 


For  a  slightly  perturbed  [0  ] , 

J 


we  then  have 
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(B-5) 


[°3]  * 


(l  +  p[NO]/[NO]o) 


1/2 


) 


which  leads  directly  to 

A[03]  (g2  -  l)  MNP] 

[°3]  "  "  2g2  [N°]o 


(  B-6) 


where  it  is  implicitly  assumed  that  ACNO]  has  a  similar  altitude  de¬ 
pendence  to  [NO]  .  Although  crude  in  its  derivation,  the  equation  above 
o 

seems  to  be  an  adequate  "rule  of  thumb."  For  g  ~  2,  it  agrees  reasonably 
with  more  detailed  model  calculations  by  Johnston23  summarized  in  Figure 
10  and  by  Chang  and  Duewer.24 
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